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Abstract 

There exists an exact relationship between the quasi-exactly solvable problems of quan- 
tum mechanics and models of square and rectangular random complex matrices. This rela- 
tionship enables one to reduce the problem of constructing topological (1/N) expansions in 
random matrix models to the problem of constructing semiclassical expansions for observ- 
ables in quasi-exactly solvable problems. Lie algebraic aspects of this relationship are also 
discussed. 



1 E-mail address: cicuta@parma.infn.it 

2 This work was partially supported by the Deutsche Forschungs Gemainschaft (DFG) and the Lodz University 
grant no. 457 

3 E-mail address: alexush@mvii.uni.lodz.pl and alexush@krysia.uni.lodz.pl 



1 



1 Introduction. The idea 



The method of orthogonal polynomials jjj is one of the basic tools for doing calculations in the 
theory of large random matrices (see e.g. recent review papers [§, ||, In order to remind 
the reader the idea of this method, let us consider for definiteness a model of random hermitean 
matrices &n of size (N + 1) x (N + 1) described by the potential V($n). The 'observables' in 
this model are initially defined as fractions of multiple integrals 

= JZWTr /($ jV )exp{-Tr V($ N )} 
{n JD$ N exp{-Tr V($ N )} U 

over components of matrix <3?at. The measure of the integration, D&n, as well as subintegral 
expressions are invariant under unitary transformations, $tv — ► Un&nU^ 1 , and this fact enables 
one to extract from both integrals (and then cancel) the volume of the unitary group by means 
of the well known Faddeev - Popov method. After this, formula (|l]) takes simpler form 

Jd\ ... d\ N A 2 (X) Elo /(A,) exp{- gjlp V{Xj)} 
IdX ...d\ N A^X)exp{-EiL V(Xi)} 

in which the variables Ao, • • • , Atv have the meaning of eigenvalues of matrix $at, and A^r(A) = 
det 1 1 A™ 1 1 is the Vandermonde determinant of these eigenvalues. 

A brilliant idea by Bessis Q was to represent this determinant in the form Ajy(X) = 
det ||P n (Aj)||, where P n (X) are polynomials normalized as P n (X) = A n + . . . and orthogonal 
with the weight e~ v ^ . This replacement enables one to reduce formula (0) to the form 

N 

(f) = E fnn (3) 
n=0 

where 

_ Jdx f{x)^ 2 n (x) 

and 

* n (x) = P n (x)e- V ^l 2 . (5) 

From the standard theory of orthogonal polynomials (see, e.g. ||) it follows that, if function 
f(x) is a polynomial of degree k, then the quantities f nn can be found from simple recurrence 
relations of the depth k + 2. In large N limit these recurrence relations can be reduced to 
a differential equation whose solution determines the leading term of the topological (1/iV) 
expansion for (|l|). The construction of the next terms of this topological expansion is generally 
rather complicated matter, because the n-dependence of numbers /„„ is not always smooth ||. 
The idea of the approach which we intend to present in this paper is to find such random matrix 
models for which the functions (|5|) could be identified with orthogonal wavefunctions for some 
quantum mechanical model. In this case the quantity (/) would have the meaning of the trace 
of function f(x) over the first N excitations in this model. In large limit, this trace could 
easily be computed (with an arbitrary accuracy) in the framework of WKB approximation. In 
this case the computation of several first terms of the topological expansion for ([!]) should not 
encounter serious difficulties. 
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Of course, such an identification is not always possible, because the wavefunctions of quantum 
mechanical models are not obliged to have the form (g). There are however, two classes of models 
for which the form of wavefunctions is given by formula (|5|). These are exactly solvable models 
associated with classical orthogonal polynomials (see e.g. and the so-called quasi- exactly 
solvable models associated with non-standard orthogonal polynomials Q. The corresponding 
random matrix models we respectively shall call exactly- and quasi-exactly solvable random 
matrix models. 

It is not difficult to see that the exactly solvable random matrix models are really exactly 
solvable in the standard sense of this word. Consider, for example, the simplest case of the 
harmonic oscillator with the potential W(x) = x 2 . The wavefunctions in this model have the 
form ^f n (x) = H n (x)e~ x / 2 , where H n (x) are Hermite polynomials. Comparing this form with 
dH) we find the form of the random matrix potential: V(x) = x 2 . But this potential describes 
the Gauss ensemble of random matrices and is exactly solvable || . It seems possible that other 
exactly solvable quantum models like Morse potential, Poschel-Teiler potential well, etc. 0, 
associated with other classical orthogonal polynomials, will be related to random matrix models 
exactly solvable for any N. 

Much more non-trivial and interesting classes of random matrix models appear if one starts 
with the quasi-exactly solvable models in quantum mechanics [0]. Remember that these models, 
which have been discovered several years ago [10, [p], [l^, 14], are distinguished by the 
fact that they can be solved exactly (analytically) only for a finite number of eigenvalues and 
corresponding eigenfunctions. It is however remarkable that the number N of exactly calculable 
energy levels is a free parameter of the hamiltonian and can be chosen arbitrarily. The aim 
of the present paper is to describe the quasi-exactly solvable random matrix models associated 
with some simplest quasi-exactly solvable quantum mechanical problems. 



2 Quasi-exactly solvable sextic oscillator model 

Consider a one-dimensional sextic anharmonic oscillator with hamiltonian 



H 



N 



h N-Q^2+ i b - [8 - 2h N ]a}x z + 4abx 4 + Aa z x b , 



(6) 



in which a > and b are real parameters and N is an arbitrarily chosen non-negative integer. 
The parameter Tin is defined as Tin = (N + 1) _1 and has the meaning of the "quantized" Planck 
constant. It is not difficult to show that the model @ is quasi-exactly solvable and has the 
order N + 1 || . This means that Schrodinger equation for (||) for any given N admits exact 
(algebraic) solutions only for N + 1 energy levels and corresponding wavefunctions. An explicit 
form of these solutions is given by the formulas 



and 



N 



E n = (A- 3h N )b + 8ah N €n,i, 



i=i 



n 



N 



(7) 



*n(a:) 



N 
1=1 



exp 



ax + bx 



2h 



N 



n 



0,...,N, 



(8) 
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in which {^i, . . . , ^ u ,n}, n = 0, . . . ,N are N + 1 different sets of real numbers satisfying the 
system of N algebraic equations 

It can be shown that N + 1 solutions of system @ describe the first iV + 1 even energy levels 
in model (@), i.e. levels with numbers 2n, n = 0, . . . ,2N. The linear span of corresponding 
wavefunctions forms a N + 1-dimensional subspace of Hilbert space which we denote hy TLn- 



3 Non-standard orthogonal polynomials 

We choose the normalization constants c n in formula (|8|) in such a way as to guarantee the 
orthonormalizability of wave functions *$> n (x): 



^> n {x)^ m {x)dx = 5 nm . (10) 

Representing wavefunctions in the form 

T / n ^ / 9n I ax 4 + bx 2 ) 

^ n (x)=V n (x 2 )exp{ (11) 



2h N 



where 

N 



i=l 

we can rewrite (flCl) as 



V n (x 2 ) =C n \{(x 2 -^) , (12) 



/+oo I ax i _|_ ^,2 1 
dxV n (x 2 )V m (x 2 ) exp < — \ = 6 nm . (13) 

From formula (|13|) it follows that V n (x 2 ) can be considered as polynomials orthogonal with the 
weight exp [—(ax 4 + bx 2 )/TiN}. Of course, V n {x 2 ) are not the classical orthogonal polynomials 
because they are of the same degree N. The general form of these polynomials is 

N 

V n (x 2 ) = J2 Pnmx 2m , n = 0, . . . , N, (14) 

m=0 

where P nm is a certain non-degenerate (and non-triangular!) (N + 1) x (N + 1) matrix. 

It is worth stressing that in the theory of random matrices a basic tool is provided by the 
set of monic polynomials (P n = x n + lower degree monomials) which are orthogonal 

/ dxP n (x 2 )P m (x 2 ) exp < - aX + bX \ = h n 5 nm . (15) 

J~oo [ h>N J 

They obey the recursion relation 

X P n {x) = P n+1 + RnPn-l (16) 
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and may be considered completely known, since Bessis yj evaluated the generating function of 
the momenta fik 

f+°° , k \ ax 4 + 

fj>k = dxx exp < — 

J-oo [ «jv 

and related the coefficients R n to the momenta. 




4 From quantum mechanics to random matrix theory 

Let us now take an arbitrary even function f{x 2 ) and consider its trace Tr^/X^ 2 ) in the space 



H.N, i-e. in the space of all exactly calculable wavefunctions: 

r+oo 



Tr N f(x' 



N 

£ 

n=0 



f(x')^ n (x)dx. 



(18) 



Using (11) and fll3|), we can rewrite (|iq ) as 



rp 2, "f^dxfix^ix^expl-h^iaxi + bx*)] 

n=0 



S-^ dxV 2 (x 2 ) expt-^^ax 4 + bx 2 )] 



(19) 



It is not difficult to see that this expression can also be represented in the form of a fraction of 
multiple integrals 



Tr N f(x 2 



■f+g IIS 1 dx % {ES 1 fjxj)} det | |Pn(x 2 )|| 2 exp {-fijg gSVl + K 2 )} 
/_ + ~ nS 1 det | \V n (x 2 )\\ 2 exp ESfM + K 2 )} 



(20) 



Using formula (14), we can write 



det N'Pnfx,- 



det 1 1 P nm Xi 



2m I 



det I IP™™ 1 1 • det I |x 



2m I 



after which formula (|2^) takes the form 



Triv/(x 2 ) 



x i "•'ft 



jV+1 

det ||P nm || • Y[ ( x i ~ x l)> 

i<k 



exp^^ESW + K 2 )} 



(21) 



1^ Ulf dx t UU 1 ^ - 4) exp Elfiaxf + bx 2 )} 

But from the random matrix theory we know that this expression can be rewritten as 

Tr N f(x 2 ) 



(22) 



J Dpi,, <t> N ] {Trf(& N $ N )} exp {-^Tr 


a($j v $ A r) 2 + b& N <S> N 


} 


jD[& N ,<S> N ]exp{-^Tr 


a(& N $ N ) 2 + b& N ® N 


} 



(23) 



where the integration is performed over all (N + 1) x (N + 1) complex matrices 3>/v It is 
well known [15] that the invariance of the subintegrals in eq.(^) under the transformation 
$ 7v — ► U\<&nU2 of the double unitary group U(N + 1) x U(N + 1) leads to the equality between 
eq.(||) and eq.(||. 
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Note also that the same steps hold if the monic classical polynomials P n replace V n in the 
above equations, then the equality of ( |20| ) with (|23| ) is a special case of evaluation of connected 
correlation functions of U(N + 1) invariant operators [16|. 



Formula fl23|) establishes an exact correspondence between the random matrix model with 
quartic potential V(& n $n) = b& N ® n + a(& N ® n) 2 (in which $ N is a {N + l) x (N + l) random 
complex matrix and $jy is its adjoint) and quasi-exactly solvable anharmonic oscillator model 
with sextic potential Wn(x 2 ) = [b 2 — (8 — 2fiM)a\x 2 + 4abx 4 + 4a 2 x 6 (in which N + 1 is the 
number of exactly calculable eigenvalues and corresponding eigenfunctions) . 

5 Rectangular random matrices 

Let us now consider another class of quasi-exactly solvable models described by the hamiltonian 

H N = -n 2 N ^ + ft|r (4S ~ 4 X ( 2 4S ~ 3) + {b 2 ~ [8 + (8s - 4)h N ]a}x 2 + 4abx A + 4a 2 x 6 (24) 

and defined on the positive half-axis x 6 [0, oo]. These models differ from the models (||) by the 
presence of an additional parameter s. If s = 1/4 then (24) reduce to (g). The exactly calculable 
wavefunctions in models (f24|) have the form 



^(x)=P n (x 2 )(x 2 r^exp|-^i-^|, (25) 

where V n (x 2 ) are polynomials orthogonal with the weight |x| 4s_1 exp[— (ax 4 + bx 2 )/%N\- Note 
however that the weight function contains now an additional factor |x| 4s_1 , and therefore, instead 
of formula (^21), we obtain 



(26) 



TW(x 2 ) = 

Jo + °° nf = r -r-'dx, {E^ 1 /(*?)} ngg(j - x 2 ) exp {-ftg gjgVI + k 2 )} 

Jo + °° nS 1 ^cfa* ngft*? - 4) ^p {-^ eS 1 W + K 2 )} 



There are two ways of transforming this expression into the fraction of matrix integrals. The first 
way is direct. We rewrite the product x t S ~ 1 as ex P ~~ 1/2) X^l 1 lnx 2 j and obtain the 

random matrix theory with a non-polynomial potential V(& n $n) = — (2s — 1/2) ln($L$jv) + 
h^ 1 [6$]y$jv + a ( < ^jv^ > A r ) 2 ]- The second way is more interesting. It is based on the observation 
that the expression ( p6[ ) naturally arises in the model of random complex rectangular matrices 
of the size (JV+l)x4s(JV+l) with quartic potential V ($%$%) = hJf 1 [b$ s N $ s N + a(® s N $ s N ) 2 ] 
(and without any singular term!). Here & S N is the transpose of the conjugate of the matrix 
This enables one to assert that relation (|23| ) holds for all quasi-exactly solvable models of the 
form p3) but the integration in the right-hand side is generally performed over all complex 



rectangular matrices of size (N + 1) x 4s (N + 1) |17], |1§, |19|, |2' 



6 The large N limit 

Up to now we considered ./V as an arbitrary non-negative integer. Now consider the case when 
./V — > oo. In large N limit the models (Q) and fl24| ) cease to be quasi-exactly solvable and become 
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exactly non-solvable. The iV-dependence of their potentials disappears and they take the form 



W(x 2 ) = {b 2 



)x 2 + 4abx 4 + 4aV. 



(27) 



The quantized Planck constant fi n tends to zero and we obtain a typical semi-classical situation 
for all the spectrum of the model (p7|). In practice, however, we do not need the information of all 
the states because the left hand side of formula (|2^) requires the knowledge of only first N + 1 
even excitations. For them we can use the standard semi-classical expressions. Substituting 
these expressions into formula (113) we obtain 



v § 



f(x 2 )dx 



Tr N f(x 2 ) = £ 



y/E 2 n-W(x 2 ) 
dx 

n=0 ^ V E 2n~W(x2) 

where Em can be found from the Bohr quantization rule 



(28) 



1 

2^ 



dx\ E- 



^2n 



W(x 2 



Hn 2n + 



(29) 



Comparing (^3|) with (p3) we arrive at the final relation 



fD[& N , <S> N ] {Trf(& N $ N )} exp {-^Tr 


a($j v $ A r) 2 + b& N $ N 


} 


fD[& N ,$ N ]exp{-±Tr 


a($]v$Ar) 2 + fe^jy^Af 


} 



N § 



f(x 2 )dx 



y/E 2n -W(x 2 ) 



dx 



(30) 



n=0 



yjE 2n -W(x 2 ) 



which can be considered as an alternative tool for doing calculations in the theory of large random 
matrices. It is worth stressing that eqs.(p9|), (|30| ) correctly describe only the leading term of 
\/N expansion for random matrix model under consideration. In order to obtain the non-leading 
terms, one should take into account the corrections to the main semi-classical approximations. 
It is interesting that, if the function f(x 2 ) is a polynomial, than the computation of several 
terms in the 1/N expansion may be performed analytically. Indeed, for a polynomial f{x 2 ), 
the main contribution to the right hand side of eq. (|30|) comes from large values of n, for which 
the analysis of the various elliptic integrals, typically appearing in semi-classical expansions, 
becomes trivial. Example of such computations will be provided in a separate paper. 



7 The general case 

It is known that, after an appropriate change of the initial variable x by a new variable t = t(x), 
the hamiltonian 



H 



N 



dx 2 



+ W N (x) 



(31) 



of any one-dimensional quasi-exactly solvable model can be represented in the form 

H N = e- Vo(t)/2 {C af3 S?S? + C a S? + C}e Voi - t)/2 , (32) 
where Sf, a = — ,0,+ are first order differential operators 
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realizing a (N + l)-dimensional representation of algebra sl(2) in the space of polynomials of 
order N. The basis of this representation is formed by monomials {l,t, . . . , t }. Formula ([3l|) 
means that, up to equivalence transformation, the hamiltonian of quasi-exactly solvable model 
is an element of the universal enveloping algebra of algebra sl(2). For this reason, the linear 
span of functions t n e- v °W 2 , n = 0,1,..., N is an invariant (N + l)-dimensional subspace 
TCn of Hilbert space and therefore the general solution of Schrodinger equation for Hm in TCn 
can be represented in the form ^> n (x) = P n (t(x))e~ v °^ t ^^ 2 , n = 0, 1, . . . ,N, where P n (t) are 
certain polynomials of degree N. Because of the orthonormalizability of wavefunctions 
the polynomials P n (t) are orthogonal with the weight e~ v ^ = (dx / dt)e~ v °^ which enables one 
to repeat the reasonings given above (see formulas (|l^) - (HH)) and derive the generalized analog 
of formula ( p3| ) 

rp f( , /^{TV/(0)}exp[-Tr V (</>)] 

in which is assumed to be a random hermitian (N + 1) x (jV + 1) matrix. Sometimes it is 
convenient to represent hermitean matrix <p in the form <J>^<1> where $ is an arbitrary complex 
matrix and is its adjoint. Then we arrive at direct generalizations of formula ([23|). 

8 Random matrix theory, virial theorems and sl(2) algebra 

Assuming that iV is finite consider function f(x) of the form 

f{x) = W N {x) + \xW' N {x). (35) 



Then, according to the well known virial theorem, the left hand side of formula ( |34|) becomes 

Tr N f(x) = Tr N H N , (36) 
where Hn is the hamiltonian of the model (^). But from (|3l| ) it immediately follows that 

Tin f{x) = C aP g aP + (N + l)C, (37) 
where is the Killing - Cartan tensor of algebra sl{2). Therefore, 
f D4>{Tx[W N i4>) +#^((/>)]}exp[-Tr V (</>)] 



fD(j> exp[-Tr V ((/>)] 



C a ^g a p + iN + l)C. (38) 



Thus we have found a purely Lie algebraic expression for random matrix integrals associated 
with quasi-exactly solvable models. 



9 Conclusion 

We have completed the exposition of our approach to random matrix models associated with 
quasi-exactly solvable problems in quantum mechanics. We see that these models are actu- 
ally simpler than general random matrix models and can be solved in a systematic way by 
constructing the semiclassical expansions for the associated quantum problems. Since the pa- 
rameter of the semiclassical expansion — the quantized Planck constant is given by the formula 
%N = 1/(N + 1), this expansion is equivalent to the topological expansion (1/N expansion) in 
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random matrix models. This fact, which we consider as the main result of the present paper, 
provides a new way for doing calculations in quasi-exactly solvable random matrix models. It 
is worth stressing in this connection that these models are not some exotic 'monsters' but are 
rather simple and ordinary looking and are rather often discussed in the literature. 

In conclusion, we would like to stress that in this paper we did not intend to present some 
systematic calculations of observables in quasi-exactly solvable random matrix models. Our aim 
was only to describe a general scheme, as to the applications, we leave them for the fortcoming 
publications (see e.g. ref. pl| ). 
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